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Abstract 


This paper deals with the calculation of the stress generated in a suspension of drops 
in the r>^®sence of a uniform electric field and a pure straining motion, te..;ing into account 
that the magnetohydrodynamic effects are dominant. It is found that the stress generated in 
the suspension depends on the direction of the ai plied electric field, the dielectric con- 
stants, the viscosity coefficients, the conductivities and the permeabilities of fluids in- 
side and outside the drois. The expression of the particle stress shows that for fluids which 
are good conductors and poor dielectrics, especially for larger drops, magnetohydrodynamic 
effects tend to reduce the dependence on the direction of the applied electric field. 


Introduction 


The study of a flow system in which the electric field and the velocity field effect 
each other is called electrohydrodynamics. The applications of electrohydrodynamics are 
numereous: cryogenic fluid management in the zero-gravity environment of space, formation and 
coelecence of solid and liquid particles, electrogasdynamic high voltage and power generation 
insulation research, physicochemical hydrodynamic, heat, mass and momentum transfer, electro- 
fluid dynamics of biological systems, and atmospheric and cloud plysics.'“^ In nature and 
industry we are quite often concerned with the properties of a fluid in which small particles 
are suspended and carried cdsout by the motion of the ambient fluid. If the average distance 
between the particles is small compared with the characteristic length of the motion of the 
suspension, one cun regard the suspension as a homogeneous fluid. The problem is to determi- 
ne the rheological properties of this equivalent fluid from the knowledge of the properties 
of the cimbient fluid. The macroscopic properties of the suspension will be referred to as 
bulk properties. When the suspension is dilute, the suspending fluid Newtonian the particles 
rigid spheres, and the particle Reynolds number sufficiently small, the suspension can be 
described in bulk as a Newtonian fluid with an effective viscosity u*= u (1 + 2.5 c) , provided 
only that the particles are not subjected to an externally applied force or couple. For the 
more general case, the non-isotropic structure of the suspension usually results in a non - 
Newtonian foirm for the bulk stress tensor.® The formulation of the problem of determining 
the stress in a suspension of particles is not straight-forward, partly due to the fact that 
the bulk stress in a suspension is not obvious. Bulk stress and other bulk properties are 
defined in terms of ensemble averages of the actual quantities. This is shown by Batchelor® 
to be equivalent to defining bulk properties in terms of volume averages provided that the 
averaging volume is chosen to contain many particles and is such that the statist. ";al proper- 
ties of the suspt. nsion are uniform over it. 

In the case of a dilute suspension, which means that the flow near any one particle is 
independent of all the other particles, the contributions to the bulk stress from the vari- 
ous particles are linearly additive. The contributions may be classified in three groups: : 
The first is a purely isotropic contribution, the second is the contribution of the deviato- 
ric stress and the third represents the contribution to the bulk stress due to the presence 
of the particles. The stress in third type of contribution is termed as "particle stress" 
and only the deviatoric part of it is significant. 

Non-Newtonian behavior of a suspension can occur in general in the following cases: ^ (i) 
Non-spherical particles can cause some directional properties.® (ii) The effect of weak 
Brownian motion.^ (iii) The effect of a couple on a particle is to rotate and to generate an 
antisymmetric part of the particle stress tensor.® (iv) The effect of the shape of a de- 
formable par-icle gives rise to non-linear stress. 'Mv) Surface tension at the boundary of 
a fluid particle or elasticity of a solid particle can cause time-dependent effects and the 
suspension exhibits visco-elastic properties. (vi) In the case of .sufficiently large 
size, the inertia forces cannot be neglected in the relative motion near one particle and 
the particle stress depends non-linearly on the bulk velocity gradient.'® 

In this paper an expression is found for the particle stress tensor of a suspension of 
drops in an electric field. It is assumed that the suspension is dilute, suspending fluid 
Newtonian, the drops spheres and the particle Reynolds number sufficiently small. The expli- 
cit form of the stress tensor depends on the electromagnetic properties of the drops and 
their surroundings > therefore, the flow due to a single drop is needed in determining the 



particle stress of a surpenslon. For this reason, consideration Is given first to f'u. flow 
due to a drop in the presence of an electric field in a pure straining motion. 

Experimentally and theoretically it has been shown by Taylor^** that a circulation can 
occuj- in a drop and its surrounding in the presence of a uniform electric field. This flow 
field set up is due to the surface charge and the tangential electric field stress over the 
surface of the drop. The flow field outside the drop given in'** is very similar to a system 
in w^;ich the field carries a uniform current. In such a case the flow field is produced by 
the rotational Lorentz force due to the distorted electric current and the associated mag- 
netic field. For fluids which are poor conductors the magn.'tic effects are negligible) and 
for fluids which are good conductors and poor dielectrics, or for larger drops the magneto- 
hydrodyamic effect may be dominant.*^ 

Following the general arguments given in^^ and the pressure and the velocity in 

the fluid outside tne drop and inside the drop which is embedded in a pure straining motion 
are determined. Since the governing equations are linear in terms of velocity and the elec- 
tric field, the effects in the cases of a drop in an electric field in the absence of a pure 
straining motion and a drop in a pure straining motion without an electric field can be su- 
perimposed. 

Since it is assumed that the suspension is dilute, the flow near one drop is independent 
of other drops and co we may use the results obtained for one droj to evaluate the integral 
in the particle stress. The expression of the particle stress has two terms ) one with an aop 
lied electric field and the other without electric field. The term which depends on electxic 
field presents a directional effect; and this makes the suspension a non-Newtonian fluid. 
This shows that the suspension of drops in an electric field cannot be represented by a vis- 
cosity that is independent of the rate-of-strain. The viscosity depends on the direction of 
the applied electric field. 

The first term of the particle stress shows the effecc in the absence of an applied elec 
trie field, and the second term denotes the additional effect due to the applied electric 
field. The second term contains two parts: one is due to the absence of the magnetic effect 
in the fluid outside the drop and the other is due to the presence of the magnetic effect. 
When the magnetohydrodynamic effect is absent, the expression of the particle stress reduces 
to that of given in.^^ 


The expression of the bulk stress 

In order to establish the relation of the bulk stress to the velocitp' and stress distri- 
butions in the fluid near individual particles, the expression for the bulk stress in the 
suspension as a volume integral is used. The average volume V is chosen to contain many par- 
ticles and is such that the statistical properties of the suspension are uniform over it. 

The hydromechanical bulk stress in the suspension is 

/o<> " • 

since the effect of inertia forces in the relative motion near a particle is neglected. The 
Maxwell bulk stress tensor may be defined in a similar manner. The bulk velocity ■ gradient is 
given by 

dv , 

V*' 

where and up are the actual values of the hydromechanical stress and velocity at any 

point X in the suspension, whether it be in the ambient fluid or inside a particle. The sur- 
face and the volume of a typ.-ical particle in V will oe denoted by A^, and respectively. 

Assuming that the ambient fluid is Newtonian with the viscositv wi the hydromechanical 
hulk stress may be written as * ^ 

* ^ Jv-ly, I ^ ^ V •'w (1) 


where the summation 

»Xi V 4-114 


is over all the particles in V, The bulk velocity gradient becomes 

iMLdV JL I r u,ni dA 

* V < > 


( 2 ) 
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where n Is an outward unit normal to Aq. Neglecting inertia effects, the following can be 
written : 


f c- j dv 

''vt 





With the addition (2) and (3), the expression (1) becomes 


(3) 




The volume V-ZV« is wholly occupied by ambient fluid and the volume Vq must be regarded as 
including the interfacial layer j and the surface will be regarded as a surface just out- 
side of the interfacial layer. The first term in tl.e expression for lit in (4) is a purely 
isotropic contribution, the second is the deviatoric stress and the third represents the 
contribution to the bulk stress due to the presence of the particles. The third term in (4) 
is called the "particle stress" and is denoted with 


It is assumed that the suspension of particles is force-free and couple-free. Since the 
exertion of a couple on the particles by external means generates an antisymmetric contribu- 
tions to che bulk stress, the bulk stress in the absence of a couple exerted on particles 
thus becomes a symmetrical tensor. 


When the ratio of the convection of charge to the conduction current, which is referred 
to as the electric Reynolds number, is much smaller than unity, the influence of the electric 
stresses on the fluid is included in the model, but there is no reciprocal effect of the moti 
on on the fields. A similar situation for the magnetic field can be considered with no effect 
of the motion on the magnetic field i and the electric current density is thus given by its 
electrostatic form. For this reason, in this paper, in view of the absence of the reciprocal 
effect of the motion on the electromagnetic field, the hydromechanical bulk stress in the 
cuspensio' is considered alone. 


Governing equations 

The magnetic induction in the fluid in and out of the drop is not negligible because of 
dynamic currents are not small enough. It is assumed that tne influence of the electric stres 
ses on the fluid is included in the model, but there is no reciprocal effect of the motion on 
the fields. Therefore, the appropriate laws of electrodynamics are essentially those of elec- 
trostatics for the electric field and the electric current density, except for the magnetic 
induction field. Under the conditions considered here, the governing ecuations of electrohy- 
drodynamics aie^^: 


e « 0 . 

(5) 

V E eO , 

(6) 

0 =CE 

(7) 

VxH.J 

(8) 

mO 

(9) 

= jmV*U ♦ 3^0'E X H 

(10) 

7-Us 0 

(11) 

■e E , is the electric field intensity, J 

the electric current density, H the magnetic 


the permeability. Throughout the paper MKS units are used. Because of the absence of fluid 
motions the termtf3fuiiH in (7) and the termd*(u«M) » M in (10) are omitted. 

The boundary conditions to be applied at the interface of a drop in an electric field are 
the following^^ 
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nx^E J ■ ® 

(12) 

n.jtfE j«o 

(13) 

n |uj -0 

(14) 

n x|uj ■ 0 

(15) 

nxf<r + t4.h} -0 

(16) 


(17) 


where CT is the viscous stress, t and h are the electric part and magnetic part of the 
Maxwell stress tensor, respectively , and denotes the jump of A across the interface. 

T is the surface tension, and Ri and R 2 are the radii of curvature of the surface (these 
radii are taken as positive when the corresponding center of curvature lies on the side of 
the interface to which n points. 

Under the conditions considered here, the electric field E and d and the velocity 
field u can be determined independently by enuatlons (5) -(11) and then, they can be related 
by the boundary conditions (12)-(17). 


A drop in an electric field 


We consider a drop or bubble, assuming that its shape is hpherlcal with radius a. The 
distance between electrodes and the drop measures to many radii thus causing the electric 
field far from the drop to be uniform. An appropriate spherical polar coordinates are defi- 
ned with the origin at the center of the drop and the symmetry axis in the direction of the 
applied field. There are four boundary conditions for the electric field intensity s (i) E 
;s finite inside the drop 1 (11) the tangential component of the electric field is continuous 
across the surface of the drop ; (lii) there is no surface current t and (iv) E tends to Bn as 
ixl tends to infinity. Subject to these boundary conditions, equations (5) -(7), (12) and (13) 
give for the electric field outside the drop 


E 




1-ot 
* 2*0t 



<-ot 


p 


(18) 


where r*.|xl , and for that inside the drop 

E = _1_ E, (19) 

2 

where . a symbol with a bar is used for the quantities of the medium inside the drop 

and a symbol without a bar is used for those of the medium outside the drop. The express.' on 
in equation (19) shows that the electric field inside the drop is uniform. 

Since there is no applied magnetic field, ecuatlons (8) and (9) give for the magnetic 
field outside the drop 

J?!)«rE, (20) 

?. 2 * ^ 7 » ' • 

and for that Inside the drop 

H =__IL_/XxE,) (21) 

The circulation in and round the drop is responsible for the electric surface force den- 
sity and the magnetic force density which are related to the Maxwell stress tensor. Expres- 
sion of t and h over the surface of the drop are needed. The tangential and normal c'>mpo- 
nent differences of t across the surface of the drop are 


n xft j - («|3-2 )(OxE.)Cii .E.) . 
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where p» f/r Is the ratio of the permittivities. The tangential and normal component diffe- 
rences of h across the surface of the drop are 


nx[h| =0 , 


n . 



9a’cr* 

* 8(24.«)* 


(X->0(n «E,)* 


The flow considered in this paper is governed by equations (10) and (11). The boundary 
conditions for the velocity are: (i) u is finite inside the drop and tends to zero as lx I 
tends to infinity > (il) u.n • 0 and u.n • 0 at the interface i (111) the tangential compo 
nent of the velocity across the drop is continuous i (Iv) tangential electric stress and 
tangential magnetic stress an'i tangential viscous stress are in balance at the Interface. 


Following the general arguments given in** and** the pressure and the velocity in the 
fluid outside the drop' * can be written as 

, j22) 

+ 9(rWb^^x,,x^xj h(r) (23) 

where ^ 


and the expressluis for R,Q,f,g and h are given in the appendix. Using the same reasoning as 
for outside the drop, the pressure and the velocity Inside the drop^* can be written as 



f? ir) ♦ 




5(r) 


(24) 


where the expression for R,Q,f,g and h are given in the appendix. 


(25) 


The balance of the normal stresses on the interface of the drop is given by equation (17). 
Since it is assumed that the interface of the drop is to be spherical the last term in equa- 
tion (17) is replaced by -2T/a. In order to find out whether the drop will become oblate or 
prolate under conditions where equation (17) is not quite satisfied, the Taylor technique*** 
is employed and it Ir assumed that a stress j applied normally to the 

surface of the drop, which )s necessary to )ceep it spherical. Replacing T in the modified 
form of equation (17) it is found that 


o a - 9£Ei(4-ft) »T i(rV E,* 

r. I. g 3 


9C-V E* 
2{2+«)* 


(3f tX) 


I 


3Xtf*rf r7 ~2m,5^_ (254. 

♦or I 



where . The expression of Tq has been given in*^ (an arlthmeti'- error in equation (30) 

in*^ is corrected, and 44 is re i laced by 14). The equilibrium geometry defends on Fq , namely 
the functional relation which is given by a,6,v,x,~,3 end a. For a detailed discussion the 
reader may' be referred to.*^ 


A drop in an electric field in a pure straining motion 

Since the governing equations are linear in terms of the velocity rnd electric field, the 
effects in the case of a droi *^n an electric fluid vithout a linear velocity at infinity an'* 
that of a drop embedded in a pure straining motion in the awsence of a uniform electric fi .. 
at infinity can be superimposed. Folloving the general arguments given in the previous pa' t 
graph the pressure and the velocity in the fluid outside the drop nan be written as 



• b|j R (r) 4- Q (r) B *»«>*>**’ 


(26) 


U< 


*i hi'*) *k*i*t 


(27) 


and Inside the drop as 

?-P. 


bji ft(r)^ bj,^ XjjX^ Q(r)♦Be^^X|^K^ 


Ui » bjj X 4 / (r ) ♦ b< j x^^ (n* b^j x^Xj K^ J5(r) b tr) « e,, jX^ Xi H (r ) 


(28) 

(29) 


where k,Q,f ,g,h,P,Q,f ,g,K,B,G,H,B,H are given in the appendix. Although a similar discussion 
to that given in the previous paragraph for the shape of the drop can be done, this is beyond 
ttie scope of this paper. 


The particle stress in a dilute suspension 

By the expression dilute suspension it is meant that the flow near one particle is inde- 
pendent of all the other particles. However, a simple modal‘s “* illustrates how surprisingly 
close the spheres are for concentrations which are numerically quite small. The relation (4) 
shows that for a dilute suspension the different particles in the volume V of suspension ma)ce 
linearly additive contributions to the particle stress and the particle stress obtained under 
these conditions is correct to the order of c (where c is the concentration of particles by 
volume). Thus, the results obtained in the previous section may be used to evaluate the third 
term in equation (4), However, some nee- ssary modifications must be made. The velocity, pres- 
sure and stress in the fluid will be written as 

. paP*p , * 

where E is a constant and ul , p' , <T,. are the disturbance quantities resulting f: jm the 
presence of the particle. The^particle stress becomes 

cpi 

"vA. 1"“ ' 1301 



since only the deviatorlc part of the particle stress is significant, the term -Fiii and si- 
milar terms are omitted. ’ 


Inserting equations (26) and (27) in equation*: (30) and using the well-)cnown Identities 

yn^nj dA . , 

Jn^n-^nifrtfdA « ) , 


in which the integration is carried over the surface of a sphere of radluc a, the following 
equation is obtained 


^<1 ( -rr-^ * 5 1? ‘ ^ ^ ‘ J 


( 31 ) 


in this equation 

»»«-«P)'-T «M,V) J 
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Since only devlatoric part of the particle stresa la algnlflcantt an laotroplc term In the 
expreaalon for the particle atreca ia not cunaide.ed. If \ goea to zero equation (31) redu> 
cea to that given in.^^ The first term of the peurticle atreaa, given by equation (31), ahuwa 
the utfect In the absence of an applied electric field, and the second term denotes the addi- 
tic.ial effect due to the applied electric field. The second term of the particle stress con- 
tains two parts t one is due to the absence of the magnetic effect in the fluid outside the 
drop and the other is due to the presence of the magnetic effect. Since it is assumed that 
the fluid outside the drop is more conductive than Inside the drop, the conductivity ratio a 
is viry much smaller than unity, and the flu.'d outside the drop is poor dielectric In compa- 
rison to that inside the drop, the dielectric -atio 8 is very much smaller than unity, then 
♦ >0 and the magnetic effect wor)(s to reduce the dei'endence of the particle stress on the 
direction of the applied electric field. 
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Appendix 

R(r).-2[f*r‘.A. , Q(r).2QL\ A . 

_J r 


A « -is* U ♦ 

iXr'a* 

|j_ 

aEjf r 

, JL (o«a. 


(4.«<)(7*S.i) 

i; 


12 ♦^«)* 

Jk (24«)*(4«-1f) »• 
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AO 
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